Mon. Not. R. Astron. Soc. 000, [|-| (2000) Printed 1 February 2008 (MN WI^ style file vl.4) 



The X-ray spectrum of a disk illuminated by ions 



H. C. Spruit^ and F. Haardt^ 

1 Max- Planck- Institut fiir Astrophysik, 

Postfach 1523, D- 85740 G arching, Germany 

^ Universita' dell'Insubria, Dipartimento di Scienze 

via Lucini 3, 22100 Corao, Italy 



Accepted xx Received ; in original form xx 



ABSTRACT 

The X-ray spectrum from a cool disk embedded in an ion supported torus is computed. 
The interaction of the hot ions with the disk increases the hard X-ray luminosity of 
the system. A surface layer of the disk is heated by the protons from the torus. The 
Comptonized spectrum produced by this layer has a shape that depends only weakly on 
the incident energy flux and the distance from the accreting compact object. It consists 
of a 'blue bump' of unComptonized soft photons and a flat high energy tail, reminiscent 
of the observed spectra. The hard tail becomes flatter as the thermalization depth in 
the cool disk is increased. Further evidence for ion illumination are the Li abundance 
in the secondaries of low mass X-ray binaries and the 450 keV lines sometimes seen 
in black-hole transient spectra. 
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1 INTRODUCTION 

The X-ray spectra F{E) of AGN as well as galactic black 
hole candidates (BHC) in their hard states are characterized 
by a power law of index ~ 1 and a high energy cutoff Ec 
around 200 keV (e.g. Mushotzky et al. 1993, Uhich et al. 
1997, Cappi et al. 1997, Zdziarski et al. 1996). Such spectra 
are well known to be describable by Comptonization in an 
electron scattering layer of optical depth r « 0.5 and tem- 
perature T ^ Ec/2. One of the classsical problems in X-ray 
astronomy is to explain why r and T should have just these 
values, with little variation between sources. Theoretical ar- 
guments can be given that Comptonization is in fact the 
most important interaction between matter and radiation 
at temperatures of 10-100 keV (e.g. Shakura and Sunyaev 
1973), for the inferred radiation energy densities near ac- 
creting black holes, but this does not tell us what the the 
thickness and temperature of the interaction region are. 

An optically thick accretion disk would produce spectra 
peaking at 1 keV and 10-100 eV for BHC and AGN, repec- 
tively. The conditions in BHC and AGN allow (at accre- 
tion rates well below Eddington) for a second form of accre- 
tion, an ion supported advection torus (Shapiro et al. 1976, 
Liang 1979, Ichimaru 19xx, Rees et al. 1982, Abaramowicz 
et al. 1988, Narayan and Yi 1994, Fabian and Rees 1995, 
Narayan et al. 1995, 1996). The ions in this flow are near 
their virial temperature, the electrons much cooler because 
of their weak interaction with the ions and their strong in- 
teraction with the radiation fleld. Such flows could produce, 
in principle, the kind of spectrum observed (Narayan et al. 



1995, 1996), but more physics must be invoked to restrict the 
optical depth and temperature of the flow to the observed 
narrow ranges (cf. Haardt 1997, Maraschi and Haardt 1997). 

1.1 Evaporating disks inside tori 

Various geometries for the accretion flow near a black hole 
have been developed, for a review see Collin (1997). One of 
the possibilities is an ion supported advection torus coex- 
isting with an optically thick accretion disk embedded in it 
(see flg. 1 in Collin). This possibility is attractive because 
the spectra of BHC in their high states show evidence of the 
simultaneous presence of an optically thick, thermal, accre- 
tion disk and a hotter component which produces a power 
law tail at higher photon energies (e.g. Mitsuda et al. 1984, 
Tanaka 1997). 

Theoretically, one would expect exchange of both mass 
an energy to take place between the disk and the advection 
torus. Heating of the disk surface by the hot ion supported 
flow above would lead to an 'evaporation' of the disk sur- 
face, feeding mass into the torus. Such evaporation has been 
studied in detail for the case of disks in Cataclysmic Variable 
systems by Meyer and Meyer-Hofmeister (1994). In the inner 
regions of the disk the mass available in the disk is smallest, 
and the energy budget potentially available for evaporation 
largest. If the mass flow from disk into torus increases with 
the energy dissipation rate, and if a steady state develops, 
one could therefore envisage a structure consisting of three 
regions: an outer one in which only a geometrically thin op- 
tically thick disk is present, inside this a composite region 
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with an evaporating disk inside a hot ion supported advec- 
tion torus, and inside this a region in which only an ion 
supported flow exists because aU disk mass has evaporated 
(Meyer and Meyer-Hofmeister 1994, Meyer-Hofmeister and 
Meyer 1999). 

Depending on details of the processes of mass and en- 
ergy exchange between disk and torus, the boundaries be- 
tween these regions may vary. It is not necessary that the 
structure is steady. The model has, in principle, sufBcient in- 
gredients to allow for variability and may perhaps be devel- 
oped further in the context of the various forms of variability 
seen in BHC. 



1.2 Energy exchange between disk and torus 

Energy exchange between disk and torus may be mediated 
by particles or by radiation. In its simplest form of interac- 
tion, the torus provides hard photons that are reprocessed 
by the disk, a model used extensively for interpretations of 
X-ray spectra (e.g. Wandel and Liang 1991, Reynolds et al. 
1994, Petrucci and Henri 1997, Gilfanov et al. 1999 and ref- 
erences therein). 

A more internally consistent model for this interac- 
tion is that of Haardt and Maraschi (1991). To simplify the 
discussion, assume that the accretion takes place predomi- 
nantly through an ion torus or a hot 'corona' above the disk 
(this assumption can easily be relaxed). The radiation pro- 
duced by the torus illuminates the disk below, which ther- 
malizes it into an approximate blackbody spectrum. These 
(soft) photons are Comptonized in the hot torus. In this 
model, approximately half the energy comes out as soft radi- 
ation and half as Comptonized photons. It correctly predicts 
the slope of the spectrum, and fixes a relation between the 
temperature and optical depth of the Comptonizing layer. 
An assumption about the temperate of of Comptonizing re- 
gion is still needed to get the right cutoff energy Ec, though 
pair ceation limits the possible temperatures. For further 
developments of this model see Haardt (1997). 

A second channel of energetic interaction is the hot pro- 
tons in the torus with temperature near the virial tempera- 
ture, Tp ~ Tv ~ 160rg/r MeV. At the distance dominating 
the energy release, r ~ 7rg, the protons thus have a tem- 
perature around 20 MeV. At this energy, they have a signif- 
icant penetration depth into the cool disk. They are slowed 
down mainly by Coulomb interactions with the ensemble of 
electrons inside their Debye sphere. The 'stopping depth', 
expressed in terms of the corresponding Thomson optical 
depth, is 



(1) 



Sttic In Alp — xip' 

(e.g. Ryter et al., 1970) where f3 = v^/c « {kTp/nipC^y^^ is 
the vertical component of the proton velocity, 9 = kT /mc(? 
is a measure of the temperature of the heated layer, InA ~ 
20 is a Coulomb logarithm, — /3^/(20), and ip the error 
function 



ip — — -= / e ^ dx 



(2) 



This formula holds for nonrelativistic temperatures; the rel- 
ativistic generalization has been given by Stepney (1983) 
and Stepney and Guilbert (1983). At low temperature kT < 



rric/mp kTp, x is small and the factor involving the error 
function can be expanded. This yields 
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2 COMPTONIZATION IN A LAYER HEATED 
BY PROTONS 

2.1 Estimating the depth of the Comptonizing 
layer 

Heating by protons yields a Comptonizing layer of thickness 
equal to the stopping depth Ts. This depth is a function 
of the temperature in the layer, by (|l]). The temperature on 
the other hand is determined by the heating and cooling pro- 
cesses, so that a consistent calculation of heating and cooling 
will yield both the temperature and the optical depth of the 
layer. With a simple estimate, we can now show that this 
will yield Ts and T in roughly the right range. 

The cooling process in the layer is the inverse Compton 
process, i.e. the energy loss electrons experience as they scat- 
ter the soft photons from the cool disk below. We assume 
that these soft photons are all (or mostly) produced by ther- 
malization of Comptonized photons from the heated layer, as 
in the model of Haardt and Maraschi (1991, hereafter HM). 
Since approximately half the Comptonized photons escape 
and the other half illuminates the thermalizing layer, the en- 
ergy flux in the soft photons at the base of the layer must be 
about the same as that in the escaping Comptonized pho- 
tons. Such a balance is possible only if the Comptonization 
is sufficiently strong. In terms of the Compton y-parameter 
y 49TS, it requires that j/ « 1. If the temperature is too 
low, the energy transfer from the electrons to the soft pho- 
tons is too low and the layer heats up until y ~ 1, and vice 
versa. Since the t/-parameter also determines the slope of 
the X-ray spectrum, the model yields a fixed spectral slope, 
which is in the range of the observed values. In HM the 
depth of the layer is a free parameter; in the present model, 
it is fixed by requiring the stopping depth Ts to be consis- 
tent with the resulting temperature T. A simple estimate is 
obtained by setting y — 1, or 8 — l/(4rs), and inserting into 
(^). This yields, assuming again that the protons are near 
their virial temperature: 

, 1/2 
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keV. 
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(5) 



(6) 



We conclude that proton illumination yields optical 
depths and temperatures in the right range, with only a 
weak dependence on the assumed distance from the black 
hole. Obviously, the estimate is rather crude, and more de- 
tailed calculations of the energy transfer from the protons 
to the electrons, as well as the Comptonization process are 
needed to test the model. In the following we make a first 
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step in this direction, by means of a radiative transfer cal- 
culation. 

2.2 Model problem 

The aim of the calculation reported below is to compute 
the temperature as a function of depth, together with the 
emergent photon spectrum from a layer heated by protons 
who deposit their energy according to (^. A preliminary 
account of the calculation has been given in Spruit (1997). 

We approximate the heating rate to be distributed uni- 
formly over a layer with depth Ts, where Ts is computed 
from the average temperature in this layer using eq. (^. I 
reality, the heating is somewhat non-uniform because of the 
dependence of the energy loss rate of the proton on both the 
proton velocity and the temperature. The approximation of 
a uniform energy input in the layer is deemed sufficient for 
the present exploratory purpose. 

In this layer, we solve the radiative transfer equation 
iteratively together with the temperature T{t), and the layer 
stopping depth Ts such that the heating is in balance with 
cooling by Comptonization of the soft photons. 

The assumptions and simplifications that go into the 
model are as follows. In the heated layer, the only photon 
process considered is electron scattering (Comptonization). 
Below this, we assume that electron scattering continues to 
be the dominant process down to some depth Tb. At Tb, the 
downward photons are assumed to be absorbed and their en- 
ergy reradiated upward as a black body spectrum. Thus, the 
gradual thermalization with depth by free-free processes is 
simplified by a step at depth Tb . The value of Ts is determined 
from the proton velocity P = v/c and the mean temperature 
in the heated layer. Obvious improvements are possible on 
these simplifications, by explicitly taking into account pho- 
ton production/destruction process, an by a more accurate 
treatment of the energy loss of the ions as they penetrate 
into the disk. 

The radiative transport part of the problem is simpli- 
fied by reducing the angular dependences to a one-stream 
model: only vertically upward and downward moving pho- 
tons are considered, and the electron distribution is similarly 
reduced from a 3-D to a one-dimensional Maxwellian dis- 
tribution. For the scattering cross section and the electron 
Maxwellian the relativistic expressions are used. The one- 
stream simplification is made for programming convenience 
only: leaving out the full angular dependences leads to very 
simple expressions. Discretization by a reasonable number 
of angles would still yield a very modest problem in terms 
of computing time, and is an obvious next step to improve 
the calculations. 

2.3 Numerical method 

The transport equation is of the form (e.g. Rybicki and 
Lightman, 1976): 

^= /d^p/d^n|^[fe(p')n(c.')x 

(1 + n{Lo)) - fo(p)n(a;)(l + n(c^'))], (7) 

where n{ijj) is the photon occupation number, /e the electron 
momentum distribution, p (respectively p') the electron mo- 
mentum, uj [ijj') the photon momentum vectors before (after) 



scattering, and $7 the scattering angle. The dependences of 
p' and uj' on (j),uJ, fl) follow from the collision kinematics. 

This equation is discretized in N^^ logarithmically 
spaced photon-energy bins. As photon energy scale we use 
ijj = hv/nicC?. As depth scale we use the Thomson optical 
depth tt. The number of depth points is fixed, but the grid 
is stretched such that at each iteration the stopping depth 
Ts is located at the same grid point. Half the grid points are 
used for the heated layer (0 < r < Tg), the other half for the 
scattering layer below. 

If and n~ are the occupation numbers of the upward 
and downward moving photons in bin i, the result is the set 
of 27V„ equations 

T = n± ^[{Bjr - B^,)nf ~ B^j] + ^ B,mf, (8) 

where Bij{T) is the scattering cross section integrated over 
the appropriate frequency and electron momenta, for scat- 
tering from bin i into bin j, in units of the Thomson cross 
section ctt. At the top boundary there are only outgoing 
photons: 

n- =0 (r = 0), (9) 

while at the lower boundary the downward photons are ther- 
malized into a blackbody spectrum of upward photons: 

nt =nBB{uJr) = l/[l-eM^/e^)] {r = r^). (10) 

Here tibb is the black-body occupation rmmber at temper- 
ature &b- This temperature follows from the condition that 
the net energy fiux at depth Tb vanishes: 

j uj^[n-{r^)-nBB{Lo,e^)] = Q. (11) 

This is a result of our assumption that the internal energy 
production rate in the cool disk can be neglected in compari- 
son with the incident proton energy fiux (generalizations are 
easily made by adding a term corresponding to the internal 
disspation in the disk). 

The condition of energy balance between the assumed 
heating rate hij) and the Compton cooling by the soft pho- 
tons is 

where the energy fiux F is given by 

F^Y.uj\nt ~n-). (13) 

i 

The equations (^) are discretized in optical depth by 
centered first order differences. The resulting set of nonlin- 
ear algebraic equations is solved by an iterative process. It 
turned out that full simultaneous linearization of the trans- 
fer equation, the boundary conditions and the energy equa- 
tion had very poor convergence properties. Instead, an iter- 
ation was done in which only the transfer equation and the 
upper boundary condition were linearized, while the lower 
boundary condition and the energy equation were dealt with 
by a modified succesive-substitution process after each itera- 
tion of the transfer equation. Convergence, however, was still 
problematic for large optical depths and for cases where the 
assumed photon energy range extended too far beyond the 
cutoff energy. For the cases reported here, where the optical 
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Figure 3. Emergent spectrum of a layer heated by protons at 
r = 7rg for three values of the energy flux 



Figure 1. Temperature as a function of optical depth in an elec- 
tron scattering layer heated by virialized protons incident with 
energy flux F (erg/cm^ /s), for r/rg = 7. Lower boundary (ther- 
malizing layer) is at t = 1. Crosses mark grid points 
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Figure 4. Dependence of the spectrum on distance from the 
compact object, for a fixed energy flux per cm^ of disk surface 



Figure 2. Variation with depth of the mean intensity J{E, t) 
of soft (0.2 keV) and hard (8keV) photons. Energy flux F = 
I.SIO^O, r/rg = 7. 

depth is not too large, on the order of 30-100 iterations were 
required for an accuracy of 10~^ in luminosity. 



3 RESULTS 

The parameters of the problem are the total energy flux F 
(per unit surface area of the disk), the velocity of the inci- 
dent protons, and the optical depth rb of the thermalizing 
lower boundary. Assuming the protons to be thermal and 
virialized, their mean vertical velocity component l3z as a 
function of the distance r from the compact object is 



The temperature as a function of depth is shown in figure ^ 
for a few values of F for r/vg = 7. These fluxes correspond 
to effective temperatures 9cft = {F / aY^'^ / {maC?) of 3 10"*, 
10"'^ and 310"'', respectively. The emergent spectrum for 
these energy fluxes is shown in figure ^ for r/rg = 7. The 
dependence of the spectrum on r/rg at a fixed F is shown 
in figure ^ The penetration depths are shown in figure ^ 
All of these cases were computed for Th = 1. 

The jump in temperature at t = Ts is a consequence of 
the fact that the heating rate jumps at Ts. The photon field is 
continuous at Ts as shown in figure ^ as expected because of 
their diffusion through the scattering layer. In order to bal- 
ance the energy input, the electrons have to loose enough 
energy by scattering soft photons. This requires them to be 
hot where the energy input rate is large, and causes the 
temperature to follow the heating rate. The jump will be 
smoothed if a more realistic model for the proton penetra- 
tion is used. If the protons are not monoenergetic and unidi- 
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Figure 5. Thomson depth of the ion-heated layer as a function 
of energy flux and distance from the compact object 



Figure 6. Spectrum integrated over the dislc (solid), compared 
with local spectrum at r = 7rg (broken). 



rectional but taken from a thermal distribution, for example, 
the transition in heating rate will be smoother. 

The similarity of the spectra, apart from shifts in ampli- 
tude and photon energy, is remarkable. The temperature of 
the heated layer increases only very weakly with increasing 
energy flux. The cutoff energy increases somewhat with dis- 
tance from the hole, but again this is dependence is rather 
weak. On account of the modest optical depth of the layer, 
the spectrum shows a prominent contribution from unscat- 
tered soft photons from the reprocessing depth rb. This peak 
is smeared out somewhat when the spectra are convolved 
over distance from the hole, at a given accretion rate. This 
is shown in figure ^. 

The final parameter of the problem is the depth of the 
thermalizing lower boundary, rb. As this depth is increased, 
the thickness of the passive scattering layer between the 
thermalizing depth and the heated surface layer increases. 
This has two effects on the spectrum. First, the strength 
of the 'blue bump' (consisting of unComptonized thermal 
photons from the lower boundary) decreases. Secondly the 
spectrum becomes flatter at the high energy end. This is 
shown in figure ^ With increasing depth of the thermal- 
izing boundary, the spectrum becomes both smoother and 
flatter. 

This can be understood as a combination of two fac- 
tors. First, the low energy photons produced at must 
travel through a scattering layer before reaching the base of 
the Comptonizing layer, and for this a gradient in photon 
density is needed. On arrival at the base of the Comptoniz- 
ing layer, the soft spectrum is therefore not a black body 
any more, but a 'diluted black body'. For a given soft fiux, 
there are fewer photons but of higher energy than if the 
thermalization were to take place directly at the base of the 
Comptonizing layer. Secondly, fewer Comptonized photons 
reach the thermalization layer, since the scattering layer in 
between effectively reflects them to some degree. The com- 
bination of these effects causes the output spectrum to be 
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Figure 7. Spectrum as a function of the depth of the thermalizing 
lower boundary Tb, for an energy flux of 2.2 lO'^* at r = 7rg 

harder than if the Comptonizing layer makes direct contact 
with the thermalizing surface. 

How much of the hardening effect remains if other radi- 
ation processes than pure scattering are taken into account 
requires more detailed calculations taking into account free- 
free and bound-free transitions. The thermalization of the 
downward traveling Comptonized photons into soft photons 
is done by these processes, and represented only crudely by 
the thermalization at a fixed depth assumed here. 



4 DISCUSSION 

With an admittedly somewhat simplified radiative transfer 
model I have shown that heating of a cool disk by protons 
from an ion supported advection torus produces X-ray spec- 
tra that are very reminiscent of the hard spectra of accreting 
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black holes. Like the Maraschi and Haardt model (and for 
the same reason), it yields approximately the right spectral 
slope, but in addition, it also predicts approximately the 
right optical depth of the Comptonizing layer and the cutoff 
energy of the spectrum. 

The temperature of the heated layer is insensitive to 
the energy flux, and stays around 40-60 keV. Instead of get- 
ting hotter at high energy flux, the incident energy is spent 
in upscattering a larger number of soft photons. As in the 
Maraschi and Haardt model, the reason for this lies in the 
energy balance condition. In order for the incident energy to 
be radiated as Comptonized flux, the Compton {/-parameter 
has to be of of the order unity, and the temperature of the 
layer adjusts accordingly. What is new in the results pre- 
sented here is that the optical depth of the Comptonizing 
layer also comes out naturally in the right range, due to the 
physics of Coulomb interaction between the virialized pro- 
tons with the electrons in the cool disk. This process also is 
fairly insensitive to the proton temperature (within the rel- 
evant range), so that the emergent spectrum is only a weak 
function of distance from the accreting object (figure ^). 

The agreement of the result with one of the most puz- 
zling features of the hard X-ray spectra of accreting compact 
objects, viz. the uniformity of the spectral shape, makes it 
likely that ion illumination plays a major role in the physics 
of these objects. 



4.1 Connection between disk and ADAF 

The ion torus (ADAF) is underluminous for its accretion 
rate, because the ions do not have time to transfer their 
energy to the radiating electrons before being swallowed by 
the hole. The presence of a cool disk illuminated by the ions 
of the torus would increase the luminosity of the system. 
Since the illumination process produces a hard spectrum, 
like the torus itself, the presence of the cool disk can still be 
compatible with the observed hard spectra, but the accreting 
system would not be as underluminous as an ADAF without 
a cool disk. If any interaction at all takes place, the ADAF 
can not be underluminous by several orders of magnitude, 
as in some proposed ADAF applications. 

The ion energy lost by illumination acts as a cooling 
agent on the hot ion torus, and limits the conditions under 
which it can exist. The hot ion supported flow has to be 
fed by a cool disk in some way or other, so that there must 
be at least a small region where the two coexist and the 
illumination process takes place. In the so-called intermedi- 
ate state in black hole accreters (e.g. Rutledge etal), a soft 
and hard component coexist in the X-ray spectrum. If ion 
tori actually exist in these systems, this is observational evi- 
dence that disk and torus can coexist with both contributing 
significantly to the luminosity. 



4.2 Lithium 

An independent observational indication for ion illumination 
is the observation of high Li abundances in the secondary 
stars of X-ray binaries (Martm et al. 1992, 1994a, 1995). The 
energy of virialized protons hitting the cool disk at r = 7rg 
(where the gravitational energy release peaks for accretion 
onto a hole) is around 50 MeV, just in the range where Li 



production by spallation of CNO elements becomes efficient. 
Since the observed secondaries are of spectral types known 
to destroy Li on a rather short time scale, a significant con- 
tinual source of Li is needed. As shown by Martm et al. 
(1994) the energetics of the accretion process is enough to 
explain the observed amount of Li on the secondary, if a 
fraction 10""^ of the Li produced in the disk finds its way 
to the secondary (in the form of a disk wind, for example). 
Another consequence of ion illumination would be Li and 
Be production by He nuclei from the virialized flow reacting 
with He nuclei in the disk. These reactions peak around 50 
Mev/nucleon, and are accompanied with emission of 7— lines 
at 431 and 478 keV. It is possible that the 7— lines observed 
sometimes around this energy (Gilfanov et al. 1991, Sun- 
yaev et al. 1992) are another signature of ion illumination 
(Martm et al. 1994a,b). 
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